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ABSTRACT 

The amount and distribution of non-gravitational energy feedback influences the global properties 
of the intra-cluster medium (ICM) and is of crucial importance in modeling/simulating clusters to 
be used as cosmological probes. AGNs are, arguably, of primary importance in injecting energy in 
the cluster cores. Comparing the observed entropy profiles, within rsoo, from the Representative 
XMM-Newton Cluster Structure Survey (REXCESS) to simulated entropy profiles without feedback, 
we estimate the amount of additional non-gravitational energy, Eicm, contained in the ICM. Adding 
the radiative losses we, then, estimate the total energy feedback, E Feedback, from the AGN's (the 
central AGN in most cases) into the clusters. The profiles for the energy deposition, AEjcm, in the 
inner regions differ for Cool-Core (CC) and Non Cool-Core (NCC) clusters; however they assume a 
similar profile after accounting for the radiation loss in CC clusters. This shows that although the 
central AGN pumps in energy, correlated with the halo mass, the amount of non-gravitational energy 
remaining in the ICM depends strongly on the amount of cooling in the central region. Turning the 
argument around we propose AEicm, the non-gravitational energy profiles for the ICM, to be a new 
and natural indicator of CC-vs-NCC clusters. The feedback energy scales with the temperature as 
EFeedback °c T^ p 57 , for the entire sample, with a scatter of « 14%. The mean non-gravitational energy 
per particle, within r^oo, is eic'M = 2.54 ± 0.74 keV. 
We use the NRAO/VLA Sky Survey (NVSS) source catalog to determine the radio luminosity, Lr, 



at 1.4 GHz of the central source(s) of our sample. For T sp > 3 keV, the Ep ee dback correlates with 
Lr, although with different normalization for CC and NCC clusters, indicating that AGNs provide a 
significant component of the feedback. Also, CC clusters show a greater radio luminosity for a given 
value of feedback energy than NCC clusters. For T sp < 3 keV, Epeedback anti-correlates with Lr] 
compared to higher temperature clusters, they are also significantly lower for the same value of Lr 
showing a lower efficiency of feedback. We then study the properties of the brightest central galaxy 
(BCG), and estimate the heating provided by them, and find a mild correlation between the BCG 
heating rate and the feedback energy. Finally, we find that mean mass deposition rate, inside the 
cooling radius, mildly correlates with the feedback energy. 
Subject headings: galaxies: clusters : general - X-rays: galaxies : clusters 
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1. INTRODUCTION 

Galaxy clusters and the intraclustcr medium (ICM) 
have long since been a topic of cosmological interest. 
Their global properties such as gas temperature, X-ray 
luminosity and SZ-flux, enable one to draw cosmolog- 
ical conclusions from s urvey s of galaxy cluste r s (e.g ., 
Reiprich fe Bohringerl (1 2003) : iVikhlinin et all (120091): 
Gladders et all (1200 7t ) : iKhedekar . Ma i umdar fe Das! 



(120101): IRozoet alJ teOlOD: ISehgal et all (|2011D : 
iBenson et all (|2011| ). However, these properties have 
been useful mostly because they can be predicted from 
models of structure formation. On the contrary, the 
detailed properties of the intracluster medium and 
their evolution with redshift are yet to be satisfactorily 
understood, since these properties are described by 
baryonic physics in addition t o the dark matter po- 
tentia l i n which it re si des (e .g. iNath fe Rovchowdhurvl 
(I2002T): IShawetalJ (f2010h: IBattaglia et al.1 (120111): 
iTrac Bode fc Ostrikenl27)Tlf ): Ichaudhuri fc Maiumdarl 
(|2011| )). Several physical processes, such as feedback 
from galaxies, including active galactic nuclei (AGNs), 
and/or radiative cooling of the ICM gas, are believed to 

anya,subha@tifr.rcs.in, biman@rri.rcs.in 



affect its X-ray properties. 

These non-gravitational processes tend to increase the 
entropy of the ICM , thereby diluting it, and conse- 
quently making it under-luminous in X-rays, especially 
in low temperature (and mass) clusters. Therefore it 
has been useful to observationally determine the entropy 
profiles of the ICM and compare with the expectations 
from various models, with or without feedback. This 
comparison of the observational entropy with the theo- 
retically predicted values allows one to estimate the de- 
gree of feedback. Following the literature, we denote 

the entropy by K = ksT '/n 2 J 3 \ where n e is the elec- 
tron number density and fcs, the Boltzmann constant. 
Note that this entropy is essentially the thermodynamic 
entropy, which we write as S, shorn of constants and 
logarithms. Entropy defined in this manner reflects the 
processes undergone by the gas, such as accretion, gas 
cooling and feedback. iVoit et ail (|2005l ) had shown that 
in the absence of any feedback and cooling processes, 
simulations tend to predict a power-law radial profile for 
the entropy outside the core, with a scaling K cx r , 
and that their AMR ( adaptive-mesh refinement ) and 
SPH ( smoothed particle hydrodynamics ) simulations 
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agreed to within ~ 10 %. In a previous paper, we have 
used the observed entropy profiles of the REXCESS sur- 
vey clusters (see below for details), and expressed the 
profile as a function not of radial distance, but of gas 
mass, since entropy per particle is a Lagrangian quantity. 
The physical processes which endow the gas with addi- 
tional entropy are also likely to move this gas around, and 
therefore it is better to view the distribution of entropy 
in gas mass rather than th e radia l distance ( iVoit et all 
IpOM k iNath fc Maiumdarl (pOlll) . After comparing the 
observed profiles with the benchmark entropy profile out- 
side the core without feedback, we determined the mean 
energy deposited per particle in the ICM, and also discov- 
ered that the profile of energy deposition has a univer- 
sal shape (|Chaudhuri. Nath fc Maiumdarl (120121 ); here- 
after referred to as CNM12). In this paper, we study 
the en tropy profi l e in t he inner as well as the outer re- 
gions. IVoit et all (|2005| ) found a large discrepancy in the 
scaled entropy profiles inside the core radii for their AMR 
and SPH simulations. While SPH simulations form clus- 
ters with almost power-law entropy distributions down to 
small radii, Eulerian simulations form much larger cores 
with the entropy distribution being truncated at signifi- 
cantly higher values. We use their SPH simulation results 
in order to determine a benchmark entropy profile with- 
out feedback inside the core, and then compare with the 
observed profiles of REXCESS clusters. 

We then discuss our results in light of the emerging sce- 
nario of feedback from radio galaxies in clusters. XMM- 
Newton and Chandra X-ray observations have shown 
that radio AGN are probably the principal agent heating 
the hot atmospheres of galaxies, clusters, and groups and 
suppressing cooling flows, reducing their strength. Su- 
pernovae are important to explain the observed metallic- 
ities in the ICM but their feedback is not strong enough 
to truncate cooling in massive galaxies. The main ev- 
idence for feedback from radio galaxies comes from the 
observations of numerous ga laxy clusters featurin g X-ray 
deficit low density cavities (| Birzan et al J 120081) . Such 
cavities have been observed to correlate strongly with ra- 
dio jets and lobes indicating that they are likely regions 
emptied of ICM by the expanding radio lobes. We find 
that the radio luminosity in the central regions is corre- 
lated with the feedback energy in the REXCESS galaxy 
clusters indicating a picture where AGN feedback plays 
a dominant role in both the CC and NCC clusters. Our 
study of the inner region of ICM also becomes important 
for the implications and connections to the properties of 
the brightest cluster galaxy (BCG). There is some degree 
of correlation between the feedback energy and the time 
averaged heating rate for the BCGS supporting a sce- 
nario of AGN feedback. We adopt a ACDM cosmology 
with H a = 70 km s" 1 Mpc~\ ft M = 0.3 and Qa = 0.7. 

2. THE CLUSTER SAMPLE 

The REXCESS survey (jBohringer et al.ll2007h is a sub- 
set of the REFLEX cluster catalog, which is a nearly 
complete flux limited cluster sa mple, covering 4.24 ste r 
in the southern extragalactic sky (jBohringer et al.ll2004j ) . 
The REXCESS sample consists of 31 local clusters ( z ^ 
0.2) , where the clusters are selected on the basis of their 
X-ray luminosity, L x = (0.407-20) x 10 44 h^ Q 2 erg s" 1 in 
the 0.1-2.4 kcV band, with no bias for any morphology 
type. This luminosity range selects clusters with a tern- 
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name (RXC) 


( keV ) 


(10 14 M Q ) 


(10 13 M Q ) 




1 


J0003. 8+0203 


3.83 


2.11 


1.99 


NCC 


2 


J0006. 0-3443 


5.24 


3.95 


4.48 


NCC 


3 


J0020. 7-2542 


5.54 


3.84 


4.06 


NCC 


1 


J0049. 4-2931 


2.87 


1.62 


1.66 


NCC 


5 


J0145. 0-5300 


5.81 


4.37 


4.85 


NCC 


6 


J0211. 4-4017 


2.08 


1.00 


.98 


NCC 


7 


J0225. 1-2928 


2.53 


0.96 


.73 


NCC 


8 


J0345. 7-4112 


2.28 


0.97 


.82 


CC 


9 


J0547.6-3152 


6.04 


4.98 


5.94 


NCC 


10 


J0605.8-3518 


4.93 


3.87 


4.63 


CC 


11 


J0616.8-4748 


4.18 


2.70 


2.86 


NCC 


12 


J0645. 4-5413 


7.23 


7.38 


10.08 


NCC 


13 


J0821. 8+0112 


2.81 


1.31 


1.16 


NCC 


14 


J0958. 3-1103 


5.95 


4.17 


4.43 


CC 


15 


J1044.5-0704 


3.58 


2.69 


3.32 


CC 


16 


J1141. 4-1216 


3.58 


2.27 


2.45 


CC 


17 


J1236. 7-3354 


2.77 


1.33 


1.21 


NCC 


18 


J1302. 8-0230 


3.48 


1.89 


1.80 


CC 


If) 


J1311. 4-0120 


8.67 


8.41 


10.69 


CC 


20 


J1516. 3+0005 


4.68 


3.28 


3.61 


NCC 


21 


J1516. 5-0056 


3.70 


2.59 


2.99 


NCC 


22 


J2014.8-2430 


5.75 


5.38 


7.19 


CC 


23 


J2023. 0-2056 


2.72 


1.21 


1.03 


NCC 


21 


J2048.1-1750 


5.06 


4.31 


5.50 


NCC 


25 


J2129.8-5048 


3.84 


2.26 


2.23 


NCC 


26 


J2149. 1-3041 


3.48 


2.25 


2.48 


CC 


27 


J2157.4-0747 


2.79 


1.29 


1.12 


NCC 


28 


J2217.7-3543 


4.63 


3.61 


4.37 


NCC 


29 


J2218.6-3853 


6.18 


4.92 


5.67 


NCC 


30 


J2234.5-3744 


7.32 


7.36 


9.87 


NCC 


31 


J2319.6-7313 


2.56 


1.56 


1.74 


CC 



TABLE 1 

Basic properties of the REXCESS clusters. The first 

COLUMN IS THE CLUSTER IDENTIFIER NUMBER, THE SECOND 
COLUMN GIVES THE REXCESS NAME OF THE CLUSTER, THE THIRD 

COLUMN IS THE SPECTROSCOPIC TEMPERATURE IN THE 
0.15 — 0.75 -R500 RANGE, THE FOURTH AND FIFTH COLUMNS SHOW 
THE TOTAL AND THE GAS MASS WITHIN i?500 ■ THE LAST COLUMN 
MARKS THE CLUSTE R AS EITHER NCC OR CC (SEE TABLE 1. IN 

IPratt et alI poiOD I. 



perature ^ 2 keV , and does not include galaxy groups. 
IPratt et all (|2010l) have noted that the REXCESS sam- 
ple is well suited to study the variation of entropy pro- 
files across a range of cluster masses, especially because 
the distances were chosen such that ^500 fell within the 
XMM-Newton field of view, which increased the precision 
of measurements at large radii. They also subdivided the 
sample into CC and NCC systems, defining the clusters 
with central density E(z)~ 2 n e _ > 4 x 10~ 2 cm -3 as CC 
systems (E(z) being the ratio of the Hubble constant at 
redshift z to its present value). Some basic properties of 
the REXCESS clusters arc shown in Table |TJ Our iden- 
tification of the cluster as CC or NC C (column 6 in th e 
Table) follows the convention used bv lPratt et al.l (|2010f) . 

In Figure [TJ we show the Dark Matter (DM) density 
profiles (Top panel) , the ICM temperature profiles (Mid- 
dle panel) and the ICM density profiles (Bottom panel) 
for a sub-sample of 4 NCC (shown in blue solid lines) 
and 4 CC clusters (red dashed lines) selected from the 
full REXCESS sample. The DM has NFW profiles with 
a concentration parameter 3.2 (for details see subsection 
13.21) . In the bottom panel, the ICM density profiles of the 
CC clusters show a central excess compared to the NCC 
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Fig. 1. — Dark matter NFW (Top panel), ICM temperature 
(Middle panel) and ICM density (Bottom panel) profiles for a sub- 
set of REXCESS clusters. The blue solid lines show NCC clusters 
while the red dashed lines show CC clusters in all the plots. The 
clusters are marked by their s pectroscopic temp erature as given 
in Table [T] (Also see Table 1 in lPratt et"aH 12010) )■ In the middle 
panel, two CC clusters (with T sp = 8.7&2.6keV are shown to have 
temperature profiles that do not fall in the inner region and have 
temperature profiles close to NCC clusters. 



clusters which leads to th e criterion, described above, 
used by iPratt et al.l (|2010f ) in marking CC clusters. In 
the middle panel, one clearly sees that NCC clusters have 
flat inner temperature profiles and 2 CC clusters have 
temperatures going down near the centre indicating cool- 
ing. However, 2 CC cluscrs (in brown dashed line) are 
seen to have flat profiles similar to NCC clusters near 
the centre, showing the simple definition of density ex- 
cess used to mark CC/NCC clusters may not be robust. 
Infact, the most natural way to identify a CC cluster is 
to identify cooling near the cluster centre which is akin 
to identifying energy loss. In a later section we propose 
a new way to identify CC clusters. 

3. ENTROPY PROFILES 
3.1. Initial entropy - radial profile 

Our goal is to compare the observed entropy profile 
(as a function of gas mass) with that expected with- 
out any non-gravitational processes affecting the ICM. 
For a baseline entropy profile, with which we will later 
compare the observed profiles, we turn to the results 
of numerical sim ulations without any radiative cooling. 
iVoit et aD (|2005| ) showed that their simulated SPH pro- 
files can be when well described, in the 0.1-1 r2oo range, 
by a median scaled profile given by the baseline power 
law relation, = 1.32(^— This mean profile has 
approximately 20% dispersion which would correspond 
to the theoretical uncertainty for our calculations. The 
corresponding relation for r 50 o is = 1.42(— — I 1,1 . 

1 ° ouu Ksaa y r 5 oo ' 

This relation was used as a benchmark entropy profile in 
CNM12 for estimating the energy deposition in the ICM 
outside the core. 

However, for the ICM inside the core, they found a 
large discrepancy in the scaled entropy profiles between 
the SPH and AMR simulations. A flat entropy core has 
been observed in the centre of non-radiative galaxy clus- 
ters in Eulerian grid codes. However this core is absent 
in lagrangian approaches such as SPH. 

Mitchell et al. 2009 have compared Eulerian simula- 
tions to the SPH simulations in the case of an idealized 
merger between galaxy clusters. They found the discrep- 
ancy between the core entrop y in the two case s, and of 
the same magnitude found bv I Voit et al.1 (120051 ) between 
their AMR and SPH codes. They find that the differ- 
ence is because of the different treatment of vortices in 
the two codes. While in the SPH case, the particles re- 
tain their initial entropy even after the merger, in the 
grid case the initial lower entropy is wiped out as a re- 
sult of the strong vortices generated during the mergers. 
Vazza (2011) found that the occurrence of a flat entropy 
core is mainly due to the hydrodynamical processes that 
are resolved in the Eulerian code, and that additional 
numerical effects influenced the entropy level to a much 
lesser degree. 

We have used the entropy profile obtained as a fit to 
the SPH simulation data points obtained by IVoit et all 
(2005) as our baseline entropy profile. Using the median 
points, we find a fit that is given by, 



ffl =0 .ll + 2.50( — ) 177 

-"■500 r 500 



(1) 



This expression is obtained as a result of fitting in the 
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Fig. 2. — The initial scaled entropy profile, K-t^/K 500, for our 
sample as function of gas mass corresponding to the scaled entropy 
profile as function of radii found in non-radiative simulations. The 
median line, shown in solid black, corre sponds to the med ian rela- 
tion obtained in the SPH simulations of lVoit et~a l. (2005) and the 
envelopes, shown in grey lines, depict the scatter about the median 
relation that they have obtained ( sec Figurc5 in Voit ct al. (2005) 
) . The vertical dashed line shows approximately the transition from 
within-core (r < 0.1r20o) to outside-core(0.1r200 < r < rsoo)- 

core region upto 0.1r 2 oo- 

There is of course an uncertainty associated with this 
relation, and we have used the (1-cr ) upper and lower en- 
velopes of the simulation results of lVoit et al.l (|2005D . in 
order to assess the implications of this uncertainty in our 
final result, and we will discuss it further below. 

3.2. Initial entropy profile with gas mass 

We first calculate the initial entropy profile (given 
above) as a function of gas mass, assuming hydro- 
static equlibirium in a Navarro-Fre nk- White (NFW) 
halo (H 

avarro, Frenk fc WhTtel fl~997l) with a concen- 
tration parameter C500 = 3.2. This value was de- 
duced for a morphologica lly relaxe d cluster sample b y 
IPointecouteau et all (120051 ) (see also lPratt et al.l (|2010D ). 
The condition for hydrostatic equilibrium states that, 



dr 



~Pg- 



GM(< r) 



± g 
A" 



3/5 



2/5 3/5 



GM{< r) 



where P a 



jksT, is the gas pressure. 



(2) 

For the bound- 



ary condition, we set the gas fraction inside the virial ra- 
dius, fg(Rvir) to Qb/Q m . We solve eqn[2]for the pressure 
profile P g , from which we determine the entropy profile 
Kth(Mg), shown in Figure [2]. The median line, shown in 
solid black, corresponds t o the median rela tion obtained 
in the SPH simulations of lVoit et al.l (|2005() given by eqn. 
[I] and the envelopes, shown in grey lines, depict the scat- 
ter about the median rela t ion th at they have obtained ( 
see Figure5 in iVoit et ahl (|lQ05i) ) . The vertical dashed 
line shows approximately the transition from within-core 
(r < 0.1r2oo) to outsidc-core(0.1r 2 oo < r < r 500 ) Note, 
that this transition happens over a range in Mg/M^ 
around the vertical line for the cluster sample. 

4. FRACTIONAL ENTROPY DEVIATION AND ENERGY 
INPUT 



We follow our previous calculation in CNM12 of en- 
ergy deposition in this paper to determine the amount 
of energy deposition associated with the entropy en- 
hancement. We recall that the amount of energy de- 
position is reflected in the quantity T h s AK/ K a h s , where 
A K = K bs — Kth, since the change in energy per unit 
mass dQ = TdS cx TAK/K, remembering the distinc- 
tion between the thermodynamic entropy and the obser- 
vational definition of entropy. 

In CNM12 it was show n that in an isochoric an d iso- 
baric processes, (see also (jLlovd-davies et aTll2000( )). 



AQ: 



AK 



(7 - l)fim p K obs 
kT obs /? 2 /3 (/? 



1) AK 



(1 - i)/xm p (/? 5 / 3 - 1) K 



obs 



(isochoric) 



(isobaric) , (3) 



The ratio of the changes in energy for a given fractional 
change AK/K Q h s and T h s , for these two cases tends to 
7 in the limit of large value of (3, and is of order ~ 1.2 
for f3 < 2. For large f3, the two estimates of energy input 
per unit mass can at the most differ by 7 ~ 1.67. Given 
the simpler expression for the case of isochoric processes, 
we use the first expression in eqn [3] for our calculation. 

We first estimate the energy per particle in the ICM, 
AE IC M{M g ) = |r o6s -A^_ ; for each cluster. Figure [3] 

shows the profiles for AEjcm /T sp , the ratio of the non- 
gravitational energy in the ICM to /c_bT sp which is in- 
dicative of the gravitational potential of the cluster. In 
CNM12 we had determined the energy deposition profile 
outside the core at Mg/M^oo > 5 x 10~ 3 , Our present 
calculations extend these to inner regions, and we show 
the individual scaled profiles, AEicm /T sp , in Figure El 
The thick blue (upper panel) and red lines (lower panel) 
show the mean profiles for NCC and CC clusters respec- 
tively. The curves show that one can clearly distinguish 
the cool core clusters from the NCC ones . The CC clus- 
ters in the inner regions dip significantly with respect to 
the NCC clusters. This is due to the energy lost by ICM 
via radiation which is estimated in later sections. 

There is however an uncertainty in the profiles that 
stems from the scatter in the benchmark entropy pro- 
files, especially inside the core, as discussed earlier. We 
show the uncertainty in Figure [4] by showing the mean 
profile, and the profiles one obtains by using the upper 
and lower envelopes of the benchmark entropy profiles. 
The mean profiles of CC and NCC clusters using the 
median theoretical relation K(r)/ K^oo(M g ) are shown 
in solid red and blue lines respectively, and the profiles 
corresponding to the lower and upper envelopes of the 
benchmark profile (shown in Figure [5]) are shown with 
dashed lines respectively. As seen in Figure [5J the me- 
dian profile, lower envelope and upper envelope differ 
only inside the core radius 0.1 r2oo- They merge out- 
side this radius. Inside the core the upper envelope in 
K t h/ is flat while the median and lower relation fol- 
low power laws down to the smallest radii. This is also 
reflected in Figure |H The upper envelope in K t h/K 500 
being much higher than the median / lower envelopes, 
the difference between the observed entropy and the thc- 



K n 



K t h reaches much 



oretical entropy i.e. AA' 
smaller values in this case. 
We then integrate the energy deposition profile and 
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Fig. 3. — Upper Panel: AEjcm /T ap profiles for all NCC clusters. 
The thick blue line is the mean line for all the NCC clusters. Thin 
lines show individual profiles for clusters with T sp < 3.5 keV (in 
green), 3.5 keV < T ap < 5keV (in magenta) and T sp > 5keV 
(in black). Lower Panel: AEicm /T sp profiles for all CC clusters. 
Individual clusters have the same colour scheme as in the top panel; 
thick red solid line shows the mean profile for all CC clusters. 

determine the total amount of energy injected, 
f kT obs AK 

ElCM= / 7 7\ ~^ dM 9> ( 4 ) 

l)(j.m p K obs 



where the integration is done from 0.05-0.5 r§Qo, to cor- 
respond to the central regions. The minimum value of 
the outermost radius upto which the data exists is 0.62 
rsoo for cluster 27. Since we have integrated all the clus- 
ter profiles over the same radial range, the upper limit 
0.5 r 50 o is chosen because all the clusters have data upto 
this radius. The lower limit 0.05rsoo is chosen since most 
cluster profiles have innermost radial points in this re- 
gion. This corresponds to typical Al g /M 500 values of 3.5 
xl0~ 4 - 0.04 . 

We show the integrated Ejcm a s a function of cluster 
temperature in Figure [7J The blue stars are for the NCC 
clusters and the red squares are for the CC clusters. The 
red dot dashed lines, blue dashed lines and the black 



Fig. 4. — The effect of theoretical (simulation) uncertainty in 
the estimation of AEjcm using EqnfT] The solid lines are for the 
median relation in K t h(r) / K^qq in simulations, with red showing 
CC and blue showing NCC clusters. The red/blue dashed lines 
show the spread in AEjcm f° r CC /NCC clusters du e to the spread 
in the theoretical entropy profiles ijVoit ct al. 2 0051 ). 

solid lines show the best fit relations for the CC clus- 
ters , the NCC clusters and the combined sample. The 
scaling relations are given in the Appendix. In CNM12 
we showed the total energy deposition Ejcm vs. cluster 
temperature, where the calculation was done for profiles 
outside the core between 0.1 r2oo (~ 0.15r5oo) and rsoo- 
Here we have looked at the energy deposition between 
0.05r5oo and 0.5rsoo- The ratio of the energy inside the 
core (i.e., r < 0.1r2oo) to the total energy within 0.5r5oo 
is ~ 9%, and to the energy within r 500 is ~ 4%; these 
are roughly the amounts by which the numbers quoted 
in CNM12 is underestimated due to the exclusion of the 
core. 

Dividing the total energy by the total number of parti- 
cles in the ICM, the mean energy per particle, e, can 
be estimated. Within rsoo, we find eicM = 2.54 ± 
0.78 keV, which is comparable to what was found in 
CNM12. If no energy had been lost in cooling, the 
ICM would have got a higher energy /particle given by 
^Feedback = 2.87 ± .75 keV . For simulations/theoretical 
modeling of AGN/SNe feedback into the ICM, the 
sources/mechanism of feedback should be able to deliver 
this energy to the gas in order to effectively model the 
energetics of the ICM. 

4.1. Input and Loss of Energy and Separation of 
Cool-Core/Non Cool- Core clusters 

In this section we take account of the energy loss due 
to cooling, which can be calculated from the observed 
X-ray bolometric luminosity, and estimate the total non- 
gravitational energy input into the ICM. This can be 
then connected to the AGN feedback which is done in 
the next section. In the process of differentiating between 
the total energy input and the non-gravitational energy 
retained in the ICM, we take another look at identifying 
CC clusters. 

4.1.1. Energetics of Cool-Core clusters 



6 



Chaudhuri, Majumdar, Nath 



Entropy 


*■ "/ 


7.3 .== 


— " > '/ 


2.8 




3 ■ 8 ____-— - 

8.7 " 












' -'2.6 






Fig. 5. — Observed entropy in CC and NCC clusters (Top panel), 
non-gravitational energy in NCC clusters (Middle panel) and non- 
gravitational energy in CC clusters (Bottom panel) for the same 
subset of REXCESS cluster as in Figure [T] Each cluster is labeled 
by their spectroscopic temperature. Top panel: blue solid lines 
show NCC clusters and red dashed lines show CC clusters. Middle 
panel: Different NCC clusters are separated by different line styles, 
with the thin linos showing AEjcm (marked by the letter 'B') and 
the thick lines showing AEp eet iback (marked by the letter 'A'). 
Bottom panel: Same as middle panel but for CC clusters. 
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Fig. 6. — The thin blue solid line shows the mean AEjcm /Tap 
profile for the NCC clusters. The thick blue solid line shows the 
mean AE Feedback /Tap profile for the NCC clusters. The red thin 
and thick lines show the corresponding profiles for CC clusters. For 
calculating AEp eedback/Tsp, tage is taken to be 2.5 Gyr. 

There have been many ways to identify cool core clus- 
ters. The REXCESS sample does it by imposing a 
minimum cutoff on the number density of electrons in 
the ICM at suitab ly defined central radii of the clusters 
(|Pratt et alJl20loh . Other ways to segregate CC/NCC 
clusters include looking at the central cooling time, the 
cuspiness in the X-ray flux profile, a ce ntral tempera- 
ture drop, value of the central entropy etc (|Hudson et al.l 
l2010f ). How ever, these diagnostic s may not be unique; 
for example iHudson et all (j2010T ) find' that "high CC" 
clusters may have a systematic central temperature drop 
whereas a "low CC" cluster may have temperature pro- 
files that are flat or only decreasing slightly towards 
the centre. This uncertainty is also seen in REXCESS 
clusters- for example, in the middle panel in Figure [1] 
two CC clusters show flat temperature profiles. However, 
the original idea of cooling flow clusters is tied to the ba- 
sic idea that the gas in the c luster cores are actually cool- 
ing w ithin a Hubble time ([Fabian. Nulsen fc Canizaresl 
119911) . Thus, for the ICM to have effectively cooled at 
any radius, the energy loss should be more than the en- 
ergy feedback. This implies that CC clusters should have 
AEjcm profiles which dip near the center and go below 
zero, implying energy taken away through cooling. Since 
AEjcm profiles directly deal with the energetics of the 
cluster, it provides a new and unique identification of 
what should be marked as a CC cluster. 

This diagnostic for cool cores is immediately seen in the 
lower panel of Figure [3] where the mean AEjc m profile 



for the CC clusters goes below zero at , , 8 ~ 10 3 . In 

° A1500 

the top panel of figure [3 we plot the observed entropy 
profiles for the same sub-sample of clusters as in figure 
Q]for comparison. CC clusters have typically lower K b s 
than NCC clusters which tend have flatter profiles at 
the center. This is due to the fact that for CC clusters 
both temperatures are lower and densities higher, thus 
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pushing the observed entropy lower than NCC clusters. 

In the middle panel and bottom panel in Figure [5l 
we show the non-gravitational energy for 4 NCC and 
4 CC clusters, respectively. The thin lines in each fig- 
ure show the AEjcm, which is the non-gravitational 
energy remaining in the ICM. As expected, all the 
NCC clusters have AEjcm positive at all radii, whereas 
for the CC clusters, AEjcm dips below zero towards 
the centre. Note, that the CC clusters with T sp = 
8.7 keV and 2.6 keV (red thin solid line and red thin dot- 
ted line), that have flat temperature profiles (see middle 
panel of Figure [1} and would be misinterpreted as be- 
ing NCC clusters if the central temperature drop is used 
as a CC diagnostic, have AEjcm dipping below zero at 
lower radii. This clearly marks these two clusters as CC 
clusters. In contrast, the CC cluster with T sp = 3.5 has 
a centrally dipping temperature profile and central den- 
sity excess but its AEjcm goes below zero at almost at 
the centre of the cluster, where the inferred ICM density 
can be influenced by the presence of the central BCG. 
Using AEicm as a NCC/CC diagnostic can also help in 
identifying CC clusters which may have been marked as 
NCC clusters. As an example, in Figurc[5J middle panel, 
we plot the AEjcm for the NCC marked cluster with 
T sp = 7.2 keV (shown in blue thin solid line with cir- 
cles) - the ICM for this cluster loses more energy than it 

recieves and has AEjcm < for ~ 10~ 3 just like 

most of the CC clusters. 
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Fig. 7. — The scaling of non-gravitational energy remaining in 
the ICM, Ejcm, after energy loss due to cooling with the cluster 
mean spectroscopic temperature T sp . Note, the relatively less en- 
ergy left for CC clusters w.r.t NCC clusters. The blue stars are for 
the NCC clusters and the red squares are for the CC clusters. The 
best-fit lines for the NCC clusters, the CC clusters and the com- 
bined sample are shown by the blue dashed line, the red dot-dashed 
line and the black solid line respectively. 



4.1.2. Energy lost due to cooling 

We estimate the bolometric luminosity Lf, i emitted by 
the ICM using the cooling function Ajv , approximated 
by Tozzi and Norman (2001) by a polynomial form for a 
mctallicity Z = 0.3Z Q . This fit reproduces the cooling 
function of Sutherland and Dopita (1993) within a few 



percent in the energy range kT > 0.03 keV. In order to 
estimate the energy lost due to cooling we need to multi- 
ply the luminosity by an appropriate time sca le, such as 
the age of the cluster. iSmith fc Tavlorl (j2008l ) have used 
semianalytic models of structure formation with lensing 
measurements to derive the ages of a sample of 10 clus- 
ters at a redshift z ~ 0.2. They derive an age of 1 Gyr to 4 
Gyr where the lower age clusters are those with disturbed 
morphology. Sinc e the redshift range of our clusters is 
similar to that in ISmith &; Tavlorl (|2008| ) , we use their 
mean age and take t age = 2.5 Gyr for our calculations. 
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Fig. 8. — The scaling of total energy feedback -^Feedback with the 
cluster mean spectroscopic temperature T sp . The blue triangles arc 
for the NCC clusters and the red circles are for the CC clusters. 
The best-fit lines for the NCC clusters, the CC clusters and the 
combined sample are shown by the blue dashed line, the red dot- 
dashed line and the black solid line respectively. 



4.1.3. Energy feedback into the ICM 

Next, we estimate the total energy deposition, 
AEF ee dback by adding the energy lost due to cooling to 
the energy which has remained in the ICM, i.e, 

AEp ee( iback = AEjcm + Lbol X t age . (5) 

Interestingly, though the mean AEjcm profiles for the 
CC and NCC clusters arc different as seen in Figure |31 
adding the energy lost due to cooling, assuming t age = 
2.5 Gyr leads to CC and NCC profiles coming within 
~ 50% of each other, the gap reducing by ~ factor of 
2 at around the core radii to more than a factor of 3 
at one-tenth the core radii. This is shown in Figure |H1 
The mean profile for AEf ee dback also falls with radius. 
Recent simulations by Gaspari et al. (2012) also show 
the energy deposition to be centrally peaked. 

For the sub-sample of 8 clusters chosen, the AEpeedback 
profiles are shown in bold lines (blue for NCC & red for 
CC, marked with the letter 'A') in the lower two panels 
in Figure [5] Notice, that once the energy lost is taken 
into account, all the profiles remain above zero. This 
input non-gravitational energy, for these clusters, finally 
becomes small close to r§oo . This is also evident in Figure 
[6l where the ratio of the feedback energy to the energy 
in the ICM due to gravity is ~ 25% at r 500 . For the CC 
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Fig. 9. — The plot shows the correlation between Bolomctric 
X-ray luminosity vs the total energy feedback for clusters in the 
REXCESS sample. The blue triangles are nore NCC clusters, with 
solid triangles for those clusters with T > 3keV and the open tri- 
angles for clusters with T < 3keV. The solid red circles are for CC 
clusters with T > 3keV and the open red circles for CC clusters 
with T < 3keV. The CC clusters with T = 8.7 keV which has 
NCC like temperature profile is shown in brown. The upper right 
quadrant, marked by the grey dashed lines contains clusters with 
T > 3 keV only, whereas the lower left quadrant has clusters with 
T < 3 keV only. 

cluster with T sp = 8.7 keV, which is one of the 8 clusters 
in the sub-sample, the energy feedback profile does not 
fall as fast as the others even outside the core. Thus, 
it has a total feedback energy which is almost an order 
of magnitude larger than that of other CC clusters of 
similar radio luminosity Lr (see, for example, Figure ITT)) . 
This cluster, typically stands out from the general trends 
seen for other CC clusters and, is separately marked in 
most of the figures. 

Finally, we can calculate the total energy Ep ee dback 
by summing over AE Feedback . This E Feedback is shown 
w.r.t. the spectroscopic temperature T sp in figure IS] The 
blue triangles are for the NCC clusters and the red circles 
are for the cool- core clusters. The red dot dashed lines, 
blue dashed lines and the black solid lines show the best 
fit relations for the CC clusters , the NCC clusters and 
the combined sample. The relations are also given in the 
appendix. The corresponding values of the scatter in the 
T- E feedback relation are 8%, 15% and 14%. The values 
for the T- Eicm relation, where Ejcm calculated using 
equation 4 does not include the cooling contribution are 
significantly higher at 29% , 15% and 22.5%. 

Figure IH1 shows the relation between E Feedback and the 
Bolometric Luminosity Lx- For a given value of Lx, 
the Epeedback values are higher for NCC clusters, with a 
typical ratio of ~ 1.7. The blue filled triangles and blue 
empty triangles are for the high temperature ( > 3keV) 
and low temperature ( < 3keV ) NCC clusters. The 
Red filled and empty circles are for the corresponding 
CC clusters. The best fit lines are shown by solid lines 
for the high temperature samples and dashed lines for 
the high temperature samples. The difference between 
the CC and NCC clusters is primarily due to the higher 
luminosity of the CC clusters for the same value of M500 / 
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Fig. 10. — This shows the radio luminosity of the central regions 
of the clusters together with their X-ray luminosity. The blue solid 
triangles are for the high temperature NCC clusters ; the blue 
hollow triangles are for the low temperature NCC clusters ; the 
red solid circles are for the high temperature CC clusters except 
for the cluster with a temperature of 8.7 keV and the brown point 
is for the CC cluster with a temperature of 8.7 keV. The blue lines 
and red lines show the best fits to the NCC clusters and the CC 
clusters. The dashed lines are for the total samples (all blue points 
or all red points plus the brown point) and solid lines are for the 
high temperature samples for both the CC and NCC clusters. 



TABLE 2 

The first column shows the cluster number 
(corresponding to cluster numbers in Table[TJ; the second, 

third and fourth columns show the radio background in 
the annulii of 15'-20', 20'-25' and 25'-30'. The fifth column 
shows the average value of the background followed by 
the background subtracted flux. note, some clusters are 
missing due to reasons mentioned in the text. 



cluster 


Llbkg 


L2 bkg 




< L hkg > 


Lr 


no. 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


1 


4.23 


4.62 


10.12 


4.45 


6.1 


2 


3.76 


3.08 


12.48 


3.50 


24.3 


3 


6.00 


2.17 


4.81 


4.33 


18.4 


6 


2.02 





0.34 


0.20 


232.1 


9 


2.96 


13.12 


2.13 


2.47 


12.7 


10 


5.78 


1.44 


3.38 


2.42 


1334.6 


13 


2.77 


3.69 


3.48 


3.34 


21.3 


15 


0.86 


2.46 


1.83 


1.78 


29.8 


16 


1.03 


0.49 


4.55 


0.80 


157.0 


18 


16.41 


2.06 


4.04 


3.01 


34.2 


19 


5.60 


4.92 


8.94 


6.49 


63.4 


20 


3.83 


4.65 


5.45 


4.64 


30.8 


21 


1.85 


1.73 


2.16 


1.91 


12.0 


22 


2.67 


2.23 


6.34 


2.51 


226.7 


26 


3.14 


1.25 


11.99 


3.77 


27.7 


27 


1.96 


8.37 


4.10 


3.02 


406.4 


28 


2.19 


3.14 


4.14 


2.67 


1.7 


29 


2.58 


4.58 


4.49 


3.79 


273.5 


30 


5.19 


1.08 


3.17 


3.20 


196.4 



T . This can be seen in the EFeedback - M500/T relations 
for the CC and NCC clusters which are very similar for 
the two sets of clusters. The typical value of the ratio 
E F eedback(NCC)/ E Fe edback(CC)(M 5 oo) ~ 1.2. Because 
of the tight relation of Epeedback with T/M500 , the low 
temperature and high temperature clusters have been 
separated by a cut on the yaxis. 
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5. ENERGY DEPOSITION AND AGN FEEDBACK 

We now compare the amount of deposited energy with 
other observed (or derived) parameters of the clusters 
that indicate the degree of AGN feedback. 

5.1. Central radio luminosity 

iMcNamara et all (|2011l) have found central radio 
sources in their sample of 400SD clusters between a red- 
shift of 0.1 and 0.6 from the NVSS 1.4 GHz catalog. 
They have considered all the sources within 250 kpc of 
the cluster center. For the background calculation they 
have used a region between 2°-4° around the center of 
the cluster. Considering the uncertainty in the centroids 
of the clusters and the NVSS resolution, they estimated 
that the sources within 250 kpc of the center are con- 
sistent with being associated with the central galaxy. In 
this manner they found radio sources within 250 kpc for 
about 30 % of the clusters. However, they did not find 
any correlation between the radio power of these sources 
and the cluster X-ray luminosity, which they found con- 
sistent with the fact that their sample is composed of 
weak and non-cooling flow clusters. The average jet 
power and the AGN heating rate that they have cal- 
culated using scaling relations was also not correlated 
with the X-ray luminosity of the clusters implying that 
the heating rate per particle was more in less massive 
systems. 

For the radio luminosity for the s ources in our sample , 
we have used the NVSS catalogue (jCondon et al.lll998l ). 
Instead of a fixed physical radius, we have taken a vari- 
able radius because the REXCESS clusters are located 
at a range of rcdshifts and the value of rsoo varies sig- 
nificantly. Therefore we have considered radio sources 
in a region 0.3 rsoo around the central BCG coordinates 
above a flux of 3mJy. We determined the background 
contribution as described below. Finally the background 
subtracted flux of the clusters in the catalog for which 
such sources exist arc listed in table [2] Clusters 5 ,8, 
11, 12, 25, 31 are not in the catalogue as they are south 
of declination —40°. Cluster 23 with 5 sources in the 
considered region is excluded. Clusters 4, 7, 17, 24 all 
of which are NCC clusters do not have a source above 
3mJy. 

We estimated the background flux by using three an- 
nulii far from the cluster center, at 15'-20', 20'-25' and 
25'-30'. We calculated the total radio flux in the annulii 
by summing the radio fluxes of all the sources in them. 
We then calculated the expected background flux in the 
region R < 0.3r5oo by scaling with the ratio of the areas. 
Table [2] gives the results of the background calculation. 
The first column is the cluster number as specified in ta- 
ble [1] The second, third & fourth columns are the total 
background (in mjy) expected in the region R < 0.3r5oo 
derived by scaling the background fluxes found in the 
annulii 15'-20', 20'-25' and 25'-30', respectively. This is 
followed by the mean value of Ll^kg, L2bk g SzL'Sbkg in mJy 
calculated using the two smaller values among the three. 
The largest value was used only if it is very close to the 
second largest. This is also because a large value in one of 
the three cases is always due a single source with a very 
high flux. The last column is the background subtracted 
flux of the radio sources within R < 0.3rsoo found by 
using this mean background value. 




Fig. 11. — The plot shows the correlation between radio lumi- 
nosity and the total energy feedback for clusters in the REXCESS 
sample. The blue triangles are for NCC clusters, with open tri- 
angles for AE calculated directly from entropy change (see text) 
. The solid triangles include the energy lost in cooling for ~ 2.5 
Gyr to show the total energy feedback. The CC clusters are simi- 
larly shown in red circles. The clusters are marked by their X-ray 
spectroscopic temperatures. The CC cluster with T = 8.7 keV 
which has a NCC like temperature profile is shown in brown. The 
clusters can be neatly separated into CC and non-core clusters oc- 
cupying two sections of plot divided by the solid black line. The 
grey dashed lines roughly group clusters in a given temperature 
range. This plot shows clusters with T > 3 keV only. 
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Fig. 12. — Same as Figure [TT] but also includes clusters with 
T < 3 keV. The different best-fit lines arc also shown for CC (T > 
3keV), NCC (both T > 3 keV and T < 3 keV) clusters (see text 
for details). Notice the opposite trend for the lower temperature 
clusters. 



We find that the background flux estimated from these 
annulii are similar except in cases where in one of the 
three annulii , the flux is significantly different because 
of one source with very high flux. In these case we use 
only the average of the other two annulii to calculate the 
background. The similarity of the background flux in at 
least two of the three annulii shows that the background 
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Fig. 13. — Plot of K-Band BCG lumi nosities with M 500 f or the 
REXCESS and the HIFLUGCS clusters H Mittal et al.ll200Sh from 
the 2MASS catalog. The REXCESS clusters are shown by the solid 
red circles and the HIFLUGCS clusters by the open blue circles. 




H (10 36 J s~ 1 ) 

Fig. 14. — The plot shows the correlation between the Heating 
Rate, H, (for details see text) vs the total energy feedback for 
clusters in the REXCESS sample. The blue triangles are for the 
NCC clusters, with solid triangles for those clusters with T > 3keV 
and the open triangles for clusters with T < 3keV. The solid red 
circles are for CC clusters with T > 3keV and the open red circles 
for CC clusters with T < 3keV. The CC cluster with T = 8.7 keV 
which has a NCC like temperature profile is shown in brown. The 
dashed line shows roughly the T = 3keV divide. 



flux converges at these radii. 

A few previous studies have found that there is a rela- 
tion between the radio flux of the central r adio source(s) 
and th e X-ray luminosity of CC clusters. I Mittal et al.l 
(|2009() found that for their strongly CC clusters, there is 
a trend between the bolomctric X-ray luminosity of their 
clusters and the radio luminosity of the central sources, 
albeit with considerable scatter. Thus there appears to 
be a connection between the AGN feedback in the clus- 
ters and the global properties such as the X-ray lumi- 
nosity. They found a correlation coefficient of 0.64 for a 
power law fit between the two quantities with the proba- 




0.7- 



L x (10' 



'Us" 1 



10 



100 



Fig. 15. — The plot shows the correlation between the bolomctric 
X-ray luminosity and the Heating Rate (for details see text) for 
clusters in the REXCESS sample. The solid blue triangles are 
for NCC clusters with T > 3kcV and the open blue triangles for 
NCC clusters with T < 3keV. The solid red circles are for CC 
clusters with T > 3keV and the open red circles for CC clusters 
with T < 3keV. The CC cluster with T = 8.7 keV which has a NCC 
like temperature profile is shown in brown. The errors on Lx are 
negligible. The blue lines and red lines show the best fits to the 
NCC clusters and the CC clusters. The dashed lines and solid 
lines are for the total samples and the high temperature samples 
for both the CC and NCC clusters. 



bility of no correlation being 2.7 x 10~ 4 . iMcNamara et all 
(2011) for a sample of 400SD clusters with a lack of 
strong CC clusters, however, did not find this correla- 
tion. 

Figure [10] plots the radio luminosity of the central re- 
gions and the X-ray luminosity of the clusters. The cor- 
relation coefficient for a power law relation between Lr 
and Lx for all the high-temperature (> 3keV) clusters 
is 0.64 with a probability of chance correlation 0.01 . 
The corresponding values for the high temperature NCC 
sample are 0.71 and 0.03 . The correlation coefficient 
for the CC sample is 0.71 with a probability of chance 
correlation 0.11 owing to the lesser no of clusters. 

Next we plot the amount of deposited energy E Feedback 
and Eicm against Lr in Figures QT] and Q21 where the 
first figure plots only the clusters with T sp > 3 keV, and 
the latter figure plots all clusters. The figures are di- 
vided with horizontal dashed lines to mark the regions of 
clusters with different ranges of T sp , as determined by the 
mean relation between AE and T sp obtained earlier. The 
CC clusters are marked red and NCC clusters are shown 
in blue. Cluster 19 with T sp ~ 8.67 keV is marked with 
a different colour as discussed before. Also we show the 
values of Ejcm and Ep ee dback (which takes into account 
the energy lost in cooling) by empty and filled symbols, 
respectively, for all clusters. The figure clearly shows that 
the CC and NCC clusters are separated in this param- 
eter space, with CC clusters occupying the right hand 
side of the plot. This is irrespective of whether one takes 
Eicm or Epeedback- In other words, for a given value of 
AE, cool clusters are more likely to have a large Lr than 
NCC clusters. 

The Figure [12] shows the same for all clusters, where 
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Fig. 16. — The plot shows the correlation between the classical 
mass deposition rate Mciasaical ( see text for details) and integrated 
total energy feedback for clusters in the REXCESS sample. The 
solid blue triangles arc for NCC clusters with T > 3keV and the 
open blue triangles for NCC clusters with T < 3keV. The solid red 
circles are for CC clusters with T > 3keV and the open red circles 
for CC clusters with T < 3keV. The CC cluster with T = 8.7 keV 
which has a NCC like temperature profile is shown in brown. The 
errors on Mciasaical are negligible. The blue lines and red lines 
show the best fits to the NCC clusters and the CC clusters. The 
dashed lines and solid lines are for the total samples and the high 
temperature samples for both the CC and NCC clusters. 
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Fig. 17. — The plot shows the correlation between the classical 
mass deposition rate M c iassical ( see text for details) and Heating 
Rate for clusters in the REXCESS sample. The solid blue triangles 
are for NCC clusters with T > 3keV and the open blue triangles 
for NCC clusters with T < 3keV. The solid red circles are for CC 
clusters with T > 3keV and the open red circles for CC clusters 
with T < 3keV. The CC cluster with T = 8.7 keV which has 
a NCC like temperature profile is shown in brown. The errors 
on M c i aas i ca i are negligible. The blue lines and red lines show the 
best fits to the NCC clusters and the CC clusters. The dashed lines 
and solid lines are for the total samples and the high temperature 
samples for both the CC and NCC clusters. 
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Fig. 18. — The plot shows the very strong correlation between the 
classical mass deposition rate M c i a ssical (see text for details) and 
Bolometric X-ray luminosity for clusters in the REXCESS sample. 
The solid blue triangles are for NCC clusters with T > 3keV and 
the open blue triangles for NCC clusters with T < 3keV. The 
solid red circles are for CC clusters with T > 3keV and the open 
red circles for CC clusters with T < 3keV. The CC clusters with 
T = 8.7 keV which has a NCC like temperature profile is shown 
in brown. The error bars along both axis are smaller than the 
marker sizes. The blue lines and red lines show the best fits to 
the NCC clusters and the CC clusters. The dashed lines and solid 
lines are for the total samples and the high temperature clusters 
respectively. 

low temperature clusters (T sp < 3 keV) arc marked 
green. The Ep ee dback for low temperature clusters (T sp < 
3 keV) is significantly lower than that of clusters above 
T sp > 3 keV for the same value of Lr showing a lower 
efficiency of feedback. We note that the trend between 
Lr and E f ee dback for the low temperature clusters runs 
opposite to that for high temperature clusters. Epeedback 
anti-correlates with Lr below this temperature i.e. clus- 
ters with a higher value of radio luminosity Lr have a 
smaller value of Epeedback or a lower efficiency of feed- 
back. However owing to the small sample size of only 
three clusters in this sample, the probability of chance 
correlation is very high. 



5.2. Correlation between BCG properties and cluster 
properties 

We turn our attention to the properties of the BCG in 
the cluster and comparisons to cluster properties. Scal- 
ing between Cluster global properties an d BCG prop- 
erties has been observed in many studies (jBrough et al.1 
10 3 20081 1 Lin fc Mohri l2004D. for which a va riety of reasons 
have been suggested. I Lin fc Mo hr (2004) have proposed 
that BCGs could grow in luminosity via mergers with 
BCGs from subclusters that have fallen into the clus- 
ter, and that the BCG would coevolve with the galaxy 
cluster. We plot the K-band BCG luminosity as deter- 
mined from the 2MASS K-band magnitudes, and M500 
of the REXCESS clusters (solid red circles) in Figure [T3l 
For comparison, we also show the da ta points (blue open 
circles) for the HIFLUGCS sample (| Mittal et al.ll2009| ) 
in the same plot . This plot shows the correlation of 
the K-band luminosity with the cluster mass M^qq. This 
naturally leads to a time averaged heating rate %{ cal- 
culated below ) that depends on the mass of the cluster 
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M500 and hence on Ef ee dback that follows a tight relation 
with M500 and T as seen in figure [H 

Next we estimate the heating rate provided by the 
BCG in the cluster and compare with the AE determined 
by us and other clu ster properties. We use the estimate 
of iBest et all ()2007l ) for the time averaged heating rate 
in terms of the stellar mass of the BCG (their eqn. 4): 

H~ 2.3 x 10 42 (A/»/10 n M Q ) ergs" 1 . (6) 

We first use the above estimated K-band luminosity of 
B CGs to determine the stellar mass, using the fits given 
bv lLonghetti fc Saraccol (|2009f ) (their eqn 9b), and finally 
obtain the heating rate using the above equation. We 
plot this BCG heating rate with AE and with L x in 
Figures [TJ] and [15] 

The correlation coefficient assuming a power law rela- 
tion between the time-averaged heating rate % and the 
Feedback energy Epeedback for the whole sample of clus- 
ters is 0.43 with a probability of chance correlation of 
0.02. The correlation coefficients for the whole CC and 
NCC samples are 0.39 and 0.45 respectively and the cor- 
responding values for the high temperature CC and NCC 
samples is 0.17 and 0.22. Thus there is a correlation be- 
tween Epeedback and the time averaged heating rate W 
when looking at the whole temperature range ( with a 
corresponding range in AE ) for CC and NCC clusters as 
well as the combined set. % is the time averaged heating 
rate due to the brightest cluster galaxies and if the BCG 
is responsible for a large fraction of the feedback, one 
would expect the energy deposited by the BCG ( oc W ) 
to be related to Ep &e dback and hence such a correlation. 

Figure 16 shows a similar degree of correlation between 
% and the bolomctric luminosity Lx ■ The correlation co- 
efficient is 0.43 with a probability of chance correlation 
0.02. This follows naturally from the above trend and 
the tight relation of the Lx and Epeedback seen in Figure 
[5] . The energy lost due to cooling and the feedback are 
expected to be of similar levels ( see for e.g. Rafferty 
et al, 2006 ). Such a relation between the the bolomet- 
ric luminosity and the heating rate is expected in this 
context. 

5.3. Correlation between the mass deposition rate and 
feedback parameters 

The classical mass deposition rate M c i ass i ca i can be de- 
rived from the density and temperature profiles of each 
cluster. It is defined as the ratio of the gas mass in- 
side the cooling radius to the cooling time t coo i i.e., 

lassical 

= M gas (r < r cool) /t coo i(r cool). The cooling- 
time tcooi is equal to the hubble time at the rcdshift of 
the cluster. Q This is a simple measure of the rate at 
which mass gets deposited and drops out of the X-ray 
band. This is true, however, if there is no source of heat- 
ing. Moreover, this should be correlated with the X-Ray 
luminosity in the core. Since, it can be argued that the 
AGN activity at the centre is influenced by the mass 
accretion rate, one expects a correlation of M c i ass i ca i 
with the radio power in the cluster centre. Note, that 
the classical mass deposition rate is not the same as the 
spectrally determined mass deposition rate, M spec which 

1 Note, the difference from the definition of fixed cooling time 
used I Mittal et al.l E009I') . 



gives the actual observed rate at which gas cools and 
feeds into the central black hole/s and should be more 
strongly correlated with the AGN output. 

6. Discussion 

There have been several studies on the dichotomy of 
CC and non CC clusters exploring the mechanisms that 
could create and destroy CCs. The two possible pro- 
cesses that have been invoked for the destruction for CCs 
are major mergers and feedback from a central AGN. 
I Burns et al.l (|2008ft have found that early major mergers 
prior to z < 0.5 are able to destroy nascent cool cores. 
Cool core clusters avoid such mergers. In their simula- 
tions which had CC and NCC clusters in the same vol- 
ume, they find that the fraction of CC clusters decreases 
with mass, as expected due to the relevance of mergers in 
the formation of higher mass clusters . However they did 
not find any dependence on redshift f Guo fc Ohl (|2009f ) 
have found that the cluster can cycle between the CC and 
NCC states with the combined effect of a time dependent 
conduction and AGN outbursts. In this scenario strong 
AGN heating can bring the CC cluster to the non CG 
state, which is maintained by conductive heating. Once 
the conduction is switched off, the cluster is cooled to 
the CC state with low level AGN heating. 

Although it appears that clusters are segregated into 
mainly two categories, CC and NCC (Cavagnolo et al 
2009), there are some cluste rs which are not easily clas- 
sified (| Rossetti et al.l |2011[ ) . The latter work found a 
connection between the presence of giant radio haloes 
and the absence of CCs. They have explored the ' evo- 
lutionary' scenario, where recent and ongoing mergers 
are responsible for the CC-NCC dichotomy. They have 
found that all clusters with radio haloes are NCC clus- 
ters, showing that processes responsible for the absence 
of the CC are also associated with the formation of radio 
haloes, thus lending support to the 'evolutionary' mech- 
anism. Combining the number of radio quiet and radio 
halo clusters with the number of CC and NCC clusters 
and assuming the lifetime of the radio halo, they find that 
the relaxation timescale of NCC clusters to CC clusters 
ranges from 1-2.7 Gyr and that this enables a NCC clus- 
ter to relax to the CC cluster. 

Rossetti & Molcndi (2010) also found regions inside 
NCC clusters that are characterized by relatively low en- 
tropy gas, and concluded that they represent the rem- 
nants of a CC after a heating event had converted the 
rest of the ICM to a NCC state. They found that in 
most cases of their sample, the heating event was related 
to merger and in a few cases, with AGN activity. 

Our results, in particular the relation between AE and 
Lfl, can be seen to support the above mentioned idea that 
CC and NCC clusters can be thought of as two different 
evolutionary stages of clusters in general. We discuss 
various aspects of this relation in light of this evolution- 
ary scenario. We recall that our main findings are: (a) 
CC and NCC clusters show a tight correlation between 
injected energy and radio (AGN) luminosity, (b) NCC 
clusters show a relatively low radio (AGN) luminosity. 

To understand these results, we point out two differ- 
ent time-scales that are relevant here: (a) the lifetime 
of r adio galaxies, wh ich is known to be tAGN < 10 s yr 
e.g. iBird et al.l 1)20081) and ( c) the cooling time for NCC 
clusters to relax to CC state, of order t coo i,NCC ~ 1 _ 3 
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Gyr (see above; I Rossetti et al.l (|2011l )). Since the first 
timescale is much shorter than the time for NCC clus- 
ters to turn to CC state, it leads to many NCC clus- 
ters without a large value of Lr corresponding to the 
amount of energy deposition. In most of these cases the 
ICM in NCC clusters have not yet reverted to the CC 
state, but the radio source is likely to have faded be- 
cause tAGN <C t C ool,NCC- This may explain the fact that 
NCC clusters in our sample have large AE but not a 
correspondingly large radio luminosity. In addition, fol- 
lowing Rossetti & Molcndi (2010), it is also possible that 
the heating is mediated by merger events, and therefore 
the NCC clusters do not show large Lr. 

7. CONCLUSIONS 

By comparing the observed entropy profiles of the 
REXCESS galaxy clusters to the baseline profile from 
nonradiative simulations we have determined the to- 
tal energy change Ejcm due to non-gravitational pro- 
cesses in line with Chaudhuri ct al 2012. The profiles 
AEicm(M 9 ) for cool core and non-cool clusters are very 
different and tend to reach negative values for cool-core 
clusters in the innermost regions due to a greater amount 
of energy lost due to radiative cooling. These profiles 
can thus be used to identify cool-core and non-cool core 
clusters in addition to other such diagnostics. Adding 
to Ejcm the energy lost due to cooling we have found 
the quantity Ep ee dback- We have studied the scaling re- 
lations of the corresponding integrated quantities Ejcm 
and EFeedback with the temperature T . These relations 



are fairly tight with the scatter in the Ejcm- T and 
EFeedback-^ relations being 22.5% and 14% respectively. 
We have calculated the radio luminosity Lr of the cen- 
tral radio sources within 0.3rsoo in the REXCESS clus- 
ters from the NVSS catalog and find that this quantity 
is correlated with the bolometric luminosity Lx for both 
CC and NCC clusters . Typically this behaviour has been 
observed in CC clusters in other studies. While the high 
temperature CC and NCC clusters show a positive rela- 
tion between Ep ee dback and Lr, the three low tempera- 
ture NCC clusters in the sample appear to show an op- 
posite trend. The Epeedback-LR relation shows a strong 
trend for both the CC and NCC clusters with similar 
power law slopes for the CC and NCC clusters. This indi- 
cates that AGN feedback from the central radio galaxies 
may provide a significant component of Epeedback hi both 
CC and NCC clusters, also supported by the correlation 
of Epeedback with the heating rate of the BCGS. The 
value of Lr in the CC clusters is however much higher 
than in NCC clusters for the same T/M 50 o/ Epeedback 
and we observe a separation in the L,R-Ep ee dback space 
for the CC and NCC clusters. 
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log 10 (E Fccdback /10 bb J) = A + Blog 10 (T/4kcV) 




A 


B 


CC 


0.13 ±0.04 


2.39 ±0.36 


NCC 


0.14 ±0.02 


2.60 ±0.15 


CC + NCC 


0.14 ±0.02 


2.56 ±0.13 


logi (Ei CM /10 55 J) = A + Blo f 


; 10 (T/4keV) 




A 


B 


CC 


0.032 ±0.055 


2.26 ±0.48 


NCC 


0.080 ±0.025 


2.58 ±0.17 


CC + NCC 


0.073 ±0.02 


2.56 ±0.16 


log 10 (E Foodback /10 55 J) = A + Blogi (L x /4 x lO^Js" 1 ) 




A 


B 


CC 


0.07 ±0.03 


0.60 ±0.08 


NCC 


0.30 ±0.02 


0.86 ±0.05 


CC + NCC 


0.22 ±0.02 


0.75 ±0.04 


CC > 3keV 


0.09 ±0.04 


0.35 ±0.11 


NCC > 3keV 


0.34 ±0.03 


0.68 ±0.09 


CC + NCC > 3keV 


0.26 ±0.02 


0.60 ±0.07 



Note: Above scaling relations are obtained using all REXCESS clusters except for cluster 14 which has much larger 
observational errors. 



loglo(EFccdback 


/10 55 J) = A + 


Blogi (L R /5 x lO^WHz- 1 ) 




A 


B 


CC 


-0.04 ±0.06 


0.22 ±0.06 


NCC > 3keV 


0.36 ±0.03 


0.22 ±0.05 


NCC < 3keV 


-0.27 ±0.01 


-0.21 ±0.12 


logi (L R /5e23WHz- 1 ) = A 


f Blogio(L x /4x lO^Js- 1 ) 




A 


B 


CC 


0.72 


0.89 


NCC 


0.38 


0.24 


CC > 3keV 


0.73 


1.36 


NCC > 3keV 


0.05 


1.69 



Note: Above scaling relations are obtained using REXCESS clusters except for (i) cluster 14 with much larger obser- 
vational errors, (ii) clusters 5,8,11,12,25,31 which are south of the region covered by the NVSS catalog, (hi) cluster 23 
which has 5 radio sources and (iv) clusters 4,7, 17 and 24 with no source over 3 mJy within .3rsoo after background 
subtraction. 



login (L x /4 x 10 37 Js- 


1 )=A- 


1- B logi (M c i assica i/50M Q yr- 1 ) 




A 


B 


CC 


-0.42 


1.09 


NCC 


-0.22 


1.13 


CC + NCC 


-0.28 


1.01 


CC > 3keV 


-0.27 


0.84 


NCC > 3keV 


.026 


0.79 


CC + NCC > 3keV 


-0.06 


0.71 
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Note: Above scaling relations are obtained using REXCESS clusters except for (i) cluster 14 and (ii) clusters 2,5, 11, 
13, 24, 25 and 27 which arc NCC clusters for which the radius r coo i is ill-defined. 



logi (E Fccdb ack/10 55 J) = A 


+ Blog 10 (H R /2x lO^Js- 1 ) 




A 


B 


cc 


0.04 ±0.03 


0.47 ±0.16 


NCC 


0.17 ±0.02 


0.63 ±0.10 


CC > 3keV 


0.08 ±0.04 


0.01 ±0.24 


NCC > 3keV 


0.17 ±0.02 


0.63 ±0.10 


log 10 (H R /2 x lO^Js- 1 ) = A + Blog 10 (L x /4 x KFjs- 1 ) 




A 


B 


CC 


-0.035 


0.055 


NCC 


0.012 


0.28 


CC > 3keV 


0.0075 


-0.10 


NCC > 3keV 


0.012 


0.25 



Note: Above scaling relations are obtained using REXCESS clusters except for (i) cluster 14 with comparatively larger 
observational errors and (ii) clusters 4 and 7 which do not have K band apparent magnitudes in the 2MASS catalog. 



log 10 (H R /2 x 10 36 Js- 


- 1 ) = A + B log^Mdassiea^OMoyr- 1 ) 




A 


B 


CC 


-.014 


-.035 


NCC 


.015 


.31 


CC > 3keV 


.052 


-.12 


NCC > 3keV 


.058 


.23 


loglo(E F ccdback/10 55 


1) = A + Blogi O (M c ia SS icai/5OM yr- 1 ) 




A 


B 


CC 


-0.10 ±0.04 


0.41 ±0.08 


NCC 


0.075 ±0.027 


0.83 ±0.082 


CC + NCC 


0.0074 ±0.02 


0.51 ±0.05 


CC > 3keV 


0.015 ±0.05 


0.21 ±0.09 


NCC > 3keV 


0.33 ±0.04 


0.51 ±0.10 


CC + NCC > 3keV 


0.21 ±0.03 


0.25 ±0.07 



Note: Above scaling relations are obtained using REXCESS clusters except for (i) cluster 14, (ii) clusters 2,5, 11, 13, 
24, 25 and 27 which are NCC clusters for which the radius R coo i is ill-defined and (ii) clusters 4, 7 for which there are 
no apparent K band magnitudes in the 2MASS catalog. 



